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Abstract 

Recent investigations have successfully 
demonstrated closed-form analytical solutions of 
spontaneous capillary flows in idealized 
cylindrical containers with interior corners. In 
this report, the theory is extended and applied to 
complex containers modeling spacecraft fuel 
tanks employing propellant management devices 
(PMDs). The specific problem investigated is 
one of spontaneous rewetting of a typical 
partially filled liquid fuel/cryogen tank with 
PMD after thrust resettling. The transients of this 
flow impact the logistics of orbital maneuvers 
and potentially tank thermal control. The general 
procedure to compute the initial condition (mean 
radius of curvature for the interface) for the 
closed-form transient flows is first outlined then 
solved for several 'complex' cylindrical tanks 
exhibiting symmetry. The utility and limitations 
of the technique as a design tool are discussed in 
a summary, which also highlights comparisons 
with NASA flight data of a model propellant tank 
with PMD. 

Introduction 

Recent investigations have successfully 
demonstrated asymptotic techniques for the 
solution of spontaneous capillary flows in 
idealized containers with interior corners 1 . The 
approach yields simple closed-form solutions for 
important features of the flow such as transient 
flow rate and 3-D interface shape without 
applying approximations such as hydraulic 
diameter, friction factors, or weighted capillary 
pressures. More recently, these techniques have 


been applied to cylindrical containers of irregular 
polygonal cross-section 2 , the results of which 
compare favorably with benchmark drop tower 
experiments. 

In this report, the theory is further extended to 
complex containers modeling spacecraft fuel 
tanks employing propellant management devices 
(PMDs). However, the general approach is 
expected to be useful to many low-gravity fluids 
management and handling operations. The 
specific problem investigated is one of 
spontaneous rewetting of a typical partially filled 
liquid fuel tank with PMD after thrust resettling. 
The transients of this flow impact the logistics of 
orbital maneuvers and potentially tank thermal 
control, particularly when the liquid inventory 
represents a significant percentage of the total 
mass of the spacecraft. 

The method of solution is briefly outlined where 
it is shown that the mean radius of curvature of 
the interface at equilibrium can be used to 
compute the pivotal initial condition for the flow 
throughout the container. This mean radius R 
may be expressed analytically for an important 
though restrictive class of simple containers 
using the approach of de Lazzer et al. 3 It is 
shown herein that this approach may be extended 
to certain more complex containers that are 
symmetric. (Computations of R using Surface 
Evolve r 4 may be employed for containers of 
arbitrary complexity.) Once R is known, the 
existing analytical solutions may be applied and 
the key characteristics of the flow may be 
determined in closed form. Examples of tanks 
with central radial and radial wall vane PMDs 
are provided. Transient flow rates are presented 
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modeling the thrust resettling problem for three 
‘complex’ containers patterned after the tank and 
PMD employed in the Vented Tank Resupply 
Experiment (VTRE Shuttle Flight Experiment). 
Despite the violation of several theoretical 
assumptions the results of comparisons to the 
VTRE data argue favorably for the use of the all- 
analytical approach as an efficient and accurate 
design tool to predict complex capillary flows in 
low-g propellant management systems. It is 
recommended that the approach also serve as a 
guide to fully transient 3-D numerical 
calculations (CFD). 

Review of Flow in ;th Comer 

Detailed comparisons between experiments and 
theory have demonstrated that spontaneous 
capillary flows in irregular polygonal containers 
with j interior corners satisfying the Concus-Finn 
corner 5 wetting condition are controlled by the 
local capillary flow in the corners 2 . Assuming a 
wetting fluid and locally parallel flow [( Hj/L ) 2 « 
1], the dimensionless leading order governing 
equations simplify to the nonlinear lubrication 
equation 

h t =2hl + hh zz , (1) 

where h = hj is the dimensionless height of the 
meniscus measured along the bisector of the yth 
corner at location z (see Fig. 1 for notation). 



Figure 1. Fluid column in an isolated comer j, angle 
2 (Xj. The 3-D surface profile is S{yzi,t) with 
characteristic height and length, H = Hj and L = Lj, 
respectively. 


This implies that the capillary surface is a 
construct of circular arcs in the cross-flow plane 
(x-y plane), and, once h(z,t) is determined, the 
entire 3-D transient surface is known from 


s j /',("/,)'(/';//• y 2 j ) x ' 2 ’ ( 2 ) 

where 

|-D | ~ h j f j sin Sj > 


and Sj =k 12-oCj-d . The parameter^ is the 

measure of interface curvature (driving force) in 
the /tli corner satisfying the Concus-Finn 
condition ( (, < n / 2 - a . ) and is given by 


fj = 


sin a j 

cos#-sinetf ; 


( 3 ) 


where 6 is the contact angle and oij is the 
particular corner half-angle. The static contact 
angle boundary condition is correct to leading 
order because the predominant flow direction is 
parallel to the contact line. The problem of 
sudden capillary rise 6 ’ 7,1 (i.e. imbibition), akin to 
termination of thruster firing during routine tank 
settling, applies constraints h(0,t) = 1, h(L,t) = 0, 
and conservation of mass to eq. (1). The solution 
for the jth interior corner provides important 
design quantities such as liquid column length Lj, 
flow rate Qj , and position of the receding bulk 

meniscus Zb as functions of time. These quantities 
are provided below in dimensional form: 


Lj =1.702 G l / 2 H) l2 t 1/2 


Qj =0.349 ffF An G j U2 H j 5l2 r y2 , 


( 4 ) 

( 5 ) 


where Hj is the constant height (a.k.a. constant 
pressure or curvature) condition at z = 0. The 
total flow rate may be determined simply as 


(L=%Qj 

j 

and the location of the receding bulk meniscus is 
approximated 2 by 


z* =i-702 


/ \l/2 

( Rcr ^ 
1 

n 


( 6 ) 


where 




-0.4103£f a „,. Fl'\ cose- since,.) 

M 


A 

R 2 




j = 1 


The geometric function 

cos 6 sin S 


J -<v 


( 7 ) 
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and Gj is given by 

g _ aF lS m 2 a j (g) 

where F- t is a weak function of 6 and a and may 
be treated as a constant P) ~ 0. 142 (see reference 
1, Fig. 6 for exact value). Note also that Fy = 
( Fi)j ~ 0. 142. a and ju are the fluid surface 
tension and dynamic viscosity, respectively. 

From eqs. (4) through (8) low gravity containers 
may be sized, fluids selected, or flow times 
predicted. Such quantities, which can be rapidly 
computed by hand, are accurate to ±6% for 
perfectly wetting fluids 6 ' 7,1 and represent an 
improvement over previous design relationships 
that used corner friction factors and weighted 
capillary pressures 8 . 

Flowever, the transients of the spontaneous 
corner flows may not be calculated without 
knowledge of H r The constant height Hj at z = 0 
is directly related to the mean radius of curvature 
R of the interface at equilibrium for the container 
in question, which is a function of container size, 
shape, fill level, and liquid contact angle(s). R 
could also depend on the fluid’s history if more 
than one local equilibrium interface 
configuration is possible. (R = U29P, where Mis 
the mean curvature of the interface.) For an 
important class of cylindrical containers with 
sufficiently planar interior corners satisfying the 
Concus-Finn condition 

R = (9) 

It is therefore necessary to determine R for the 
container before Hj and the subsequent transient 
flows in each corner may be computed. 

Calculation of Tank Mean Radius of Curvature 

In the zero-gravity environment, for cylindrical 
containers of arbitrary cross-section that possess 
at least one interior corner satisfying the Concus- 
Finn condition, de Lazzer et al. 3 apply the 
divergence theorem to the Young-Laplace-Gauss 
equation 



over a presumed solution domain O’ bounded in 
part by circular arcs of radius R that cut off 
corner flow sections and meet the rigid walls in 


the prescribed contact angle 6. See Figure 2 for 
the case of a rhombus. When such a domain can 
be found, the divergence theorem yields 

^cosff* ds = | q *|//{ ' ( 10 ) 

Zj u £ 2 

Here X! denotes the totality of boundary arcs on 
the rigid part of the boundary, and X 2 denotes the 
circular arcs that appear, see Figure 2b. On X 1; 
6 * = 0 is the prescribed contact angle of the 
liquid with the material of the containing vessel; 
in accordance with the method, 0* is set equal to 
zero on Z 2 , corresponding to the hypothesis that 
the fluid rises vertically on the X 2 arcs. 

For the cases of regular polygonal and rhombic 
cylinders, de Lazzer et al found that by inserting 
arcs symmetrically into corners as indicated in 
Figure 2, a unique value of R consistent with the 
construction could be found. We outline that 
procedure for the rhombic case in eqs. (11) to 
(13) below. It does not follow directly from the 
method that the value thus determined actually 
corresponds to a solution of the form desired; 
however the correctness of the procedure for the 
case of a regular polygon was later demonstrated 
by Finn and Neel" 1 . These authors go on to point 
out that in a general configuration the application 
of the method becomes difficult and additionally 
can lead to erroneous results. Nevertheless, the 
procedure does lead to formally solvable closed 
form expressions for R for a variety of relevant 
container section types, several of which have 
been verified experimentally: squares 1 ' 6 , rhombi 7 , 
rectangles 1 , equilateral triangles 1 , irregular 
triangles 2 , and simple cylinders with regular 
vanes 9 . Although the hazards pointed out by Finn 
and Neel are real, one may presume on the basis 
of their success with the regular polygon that at 
least some of these special cases correctly 
represent reality. Beyond that, the close 
correlation we have found in the cases we 
consider, with numerical results from the Surface 
Evolver and comparison with experiment, speak 
strongly for the underlying correctness of the 
present application. 

In a general case and especially for asymmetric 
configurations, strong caution must be advised. 
In the present paper, symmetric interfaces in 
symmetric containers will be assumed in like 
manner as in de Lazzer et al, since such 
interfaces are frequently observed in practice. 
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As illustrated in Fig. 2, for a cylinder with 
rhombic section where the Concus-Finn 
condition is only satisfied in the corners with 
acute angles 3 , along the /th portion of the 
perimeter the contact angle for use in eq. (10) is 
ft. The area contained within the projected 
perimeter is FI and is identified by a heavier line 
weight in this and figures to follow. The sought 
mean radius of curvature of the interface is R. 
The dashed lines sketched in Fig. 2b will be 
discussed shortly. 



Figure 2. Rhombic cylinder with wetting of acute 
edges only, after de Lazzer et al 3 . 


A j =R 2 F Anj =f? h }F Anj . 

Note that E of eq. (13) bears no relation to E,- of 
eq. (11). 

Modified Approach to Calculate R 

An alternative application of the technique of de 
Lazzer et al may be pursued by identifying and 
analyzing symmetric sub-sections of a given 
container cross-section. For example, the 
smallest symmetric subsection of the rhombic 
cylinder example of Fig. 2 is the quarter section 
identified by dashed lines in Fig. 2b. This 
symmetric subsection is redrawn in Fig. 3. An 
additional angle ft must be specified along the 
symmetry boundaries. Assuming the Concus- 
Finn condition is satisfied only at the acute 
vertex, eq. (11) for the geometry of Fig. 3 
becomes 

|E,|cosft + |E 2 |cosft, + |E 3 |cosft = h-|/K . (14) 


The left hand side of eq. (10) may be evaluated 
and represented as the summation of projected, 
interface perimeter lengths E,- , weighted by 
cos 6; , and enclosing area Q : 

n 

^E ; cos#, = ■ (11) 

i=l 

For the polygonal section depicted in Fig. 2, ft = 
6 ; following de Lazzer et al ft = 0, and eq. (11) 
becomes 

|E 1 |cose + |E 2 | = b*|//J, (12) 


which when solved for R yields 
PcosO 

R = - 


2 E 



f \ ,/2 l 


4AE 


1- 

1 



D 2 2 a 

P cos O 


_ 

V z 

_ 


(13) 


where P and A are the total perimeter and area of 
the container cross-section, respectively, and 


^ =5 X, 

j=i 

with F An j given by eq. (6). For the rhombic 
section of Fig. 2 is E = 2 F An . F Anj is the 
dimensionless geometric constant of proportion- 
ality for the cross-flow area Aj and mean radius 
of curvature squared; namely. 


Along the exposed (unwetted) faces of the 
rhombus 2j, ft = ft the contact angle of the 
liquid on the wall material. Along the fluid 
interface spanning the corner 2), ft 2 = 0. 
Additionally, because the dashed lines identify 
planes of symmetry for the surface, along 2), 
ft = JiH. Substitution of these quantities into 
eq. (14) produces 

|Ej|cosft +|E 2 | = |d*|//? , (15) 

which is identical to eq. (12) only 2j in this case 
does not include the symmetry plane portions of 
the perimeter of the subsection. Solving eq. (15) 
for R in this case yields 


R = 


P w cos 0 
2 E 



/ _ \" 2 1 


4AE 


1- 

1 ; — 



d 2 2 a 

P cos & 


_ 

\ w J 

_ 


(16) 


which produces the same value for R as 
computed by eq. (13) since for this symmetric 
subsection E = F An /2, and P w and A are 25% the 
values for the full domain solution, eq. (13). P w is 
the perimeter of the section minus the symmetry 
boundaries. 
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Fig. 3. Symmetric subsection for rhombus of Fig. 2. 
Symmetry planes identified by dashed fines. 

As will be demonstrated, this modified approach 
to compute R is useful in determining flows in 
more complex containers. But the technique is 
fundamentally limited by the assumption of 
symmetric interfaces in symmetric containers. 
Uniqueness and stability of particular presumed 
interfacial configurations based on intuition and 
experience may also be difficult to establish 10 
and will depend on fluid fill level and history for 
real systems. 


Calculation of R in Complex Cylindrical 
Containers with Symmetry 

Cylindrical Tank with Central Radial Vanes 


By viewing more complex container cross- 
sections as collections of symmetric subsections 
it is possible to compute R analytically for 
a variety of important container types 
with applications to low-g propellant/cryogen 
management. 

For example, a cross-section of a long, partially- 
filled, right circular cylindrical propellant tank 
model with central radial vane structure is 
sketched in Fig. 4a. Again, due to the symmetry 
of the tank the equilibrium mean radius of 
curvature of the interface R may be determined 
by analyzing the smallest symmetrical element of 
the section as sketched in Fig. 4b. Assuming the 
Concus-Finn condition is satisfied between each 
of the vanes, eq. (16) for the geometry of Fig. 4b 
yields again 


nondimensional quantities !R = R/r and 
1 / = V tr , eq. (17) becomes 


%.= 


(• V + a) 



/ r, \ 

1 / 2 “ 

1- 

y 01 F An 



y (V + CC) 2 ^ 



■ ( 18 ) 


Eq. (18) is constrained by at least the condition 
H < “Esina /sin <5 ; the interface cannot pin on 
the vane edges. Other constraints are possible, 
such as the case of wetting between the vanes 
and the circular tank wall which is not considered 
here though increasingly likely as the vane length 
^approaches 1. 

The symmetrical tank sketched in Fig. 4 may be 
generalized to a tank possessing «-vanes. For 
such a tank, and for 6=0, eq. (18) is presented 
in Fig. 5 for a variety of dimensionless vane 
lengths U. The domain of each curve is limited 
by the constraint of no pinning on the vane 
edges. As is observed from the figure, the case of 
only 2 vanes with V = 0 recovers the correct 
solution of the right circular cylinder without 
vanes, r J{= 0.5. It is also observed from the figure 
how Hi decreases with increasing number of 
vanes (decreasing a). 




Figure 4. Simplified cylindrical tank model with 
central radial vane PMD: a. cross-section identifying 
wetted vanes, b. symmetric element of shaded region 
in a. with X 3 identifying symmetry planes. 


R = 


P w cos 6 
2 X 



/ „ \ 

1/2 


4 A X 


1 - 

1 


d 2 2 a 

r, cos 6 


_ 

\ w V 

_ 


, (17) 


where I\, = 2V + 2 ar, A = a r, and X = as 
given by eq. (6) for the wetted corner formed by 
the vane of vertex angle 2a. For 6=0 , defining 


Despite the limitation of no pinning allowed on 
the vane edge, the dimensional mean radius of 
curvature of the interface R = r% may be 
computed from eq. (18) for a number of vane 
lengths V = rV of practical importance. 
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Figure 5. t^vs. “Ffor «-vaned tank patterned after the 
tank of Fig. 4. 


Regular n - gon Tanks with Radial Wall Vanes 

Another benchmark tank model readily 
addressed by the analysis outlined herein is that 
of regular polygonal cylindrical tanks with radial 
vanes emanating from the corner vertices. 
Several such tanks are sketched in Fig. 6 for 
n = 3, 4, 6, and 12. As n increases this tank 
model approaches that of a right circular 
cylindrical tank with radial vanes emanating from 
the tank wall. The tank with n = 12 is presented 
in Fig. 7 in greater detail. 




(V + 


2 F. 


1 - 


1 - 


F An sin(2 tt/ti) 
(' V + sin (Jc/n)Y 


xl/2 


(19) 



Figure 7. Regular polygonal tank with radial wall 
vanes, n = 12: a. cross-section identifying wetted 
vanes, b. symmetric element of shaded region in a. 
with Ej identifying symmetry planes. 


Eq. (19) is constrained by at least 2 conditions: 

1. % < 'Fsin a /sin 6 , interfaces can not pin on 
vane edges. 

2. % < sin(;r In) sin a / sin S , a single inter- 
face can not span two corners. 



Figure 6. Regular polygonal tanks with radial wall 
vanes: n = 3, 4, 6, and 12. 

Again, due to the symmetry of the tank the mean 
radius of curvature of the interface R may be 
determined by analyzing the smallest 
symmetrical element of the section as sketched in 
Fig. 7b for the case n = 12. This element is a 
right triangle with acute vertex angles Jt/n and 
4T/2 - 1 In). Assuming 8=0, the Concus-Finn 
condition is satisfied in each interior corner 
formed by the vanes, and eq. (17) for this 
problem may be solved for R and 
nondimensionalized by tank circumscribing 
radius r yielding 


Again, other constraints are possible, such as the 
case of a single interface wetting two vanes near 
the tank axis for large V. This case is not 
considered here though increasingly likely as the 
vane length V approaches 1. 



Figure 8. l^vs. ‘f'from eq. (19) for tank of Figs. 6 and 
7 with n = 3, 4, 6, 8, 12, 24, and 48. Dashed line at 
curve terminus implies vane pinning constraint, dot- 
dashed hne implies single interface spans two corners. 


t^is computed via eq. (19) as a function of 'Ffor 
a variety of n and presented in Fig. 8. The 
domain of the solutions is limited by at least the 
two constraints identified on the figure. For the 
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case n = 12, the curve identifying the complete 
range of with noted constraints is presented 
for later discussion. 


nQm 

nV 


0.349- 


V 


oF, R sin" a 
Mf 2 t 


\l/2 


( 22 ) 


Solution to Transient Flows using 7? 


Once 3(,is known for the tank, Hj values for each 
corner flow are computed using eq. (9) and the 
design quantities provided in eqs. (4), (5), and 
(6) may be determined. In addition, the entire 
surface profile of the liquid throughout the 
container may be computed. The solution follows 
from a global similarity solution and is 
applicable at long times throughout the container, 
despite the fact that both the flow and interface 
shape are not known in the neighborhood of the 
bulk meniscus 2 . By approximating the global 
similarity solution for the meniscus centerline 
height in each corner by the polynomial 

hj = Hj(i-0.51\lfj -0.429/jf 2 ) (20) 


with 


7/J =0.587 


Mfj 

° Hj (F/) - sin 2 (Xj 




zr 1/2 (21) 


subject to the constraint 


n b fj 
Fj 1 sin a j 


</?;<i 


the 3-D transient interface in each corner may be 
computed via eq. (2). 

Note that for the tanks of Figs. 4 and 7, the index 
j is somewhat superfluous since all interior 
corners of the tank are identical. 


Examples of Design Utility 

Cylindrical tanks may be designed with optimal 
characteristics using the analytical solution 
approach. A hypothetical example might be a 
PMD which would minimize tank rewetting time 
following resettling without an excessive mass 
penalty for unnecessary vanes. To address this 
optimization problem one might compute a ratio 
of total flow rate to total vane length. For the 
specific case of the central radial vane tank 
model sketched in Fig. 4 this ratio employs eq. 
(5) and is given by 


where V is the vane length and F A „, R, a, and / 
are functions of the number of vanes n. 
Substituting % from eq. (18) into (22) and 
retaining only dimensionless geometrically- 
dependent terms, one computes 


0,=- 


100 sin« 
fFUFV 


f 

r r 

jCF+a) 

i 

l 


1— 


a Fa,. 
(F + a ) 2 


,1,2 


7 JJ 


(23) 


where the prefactor of 100 serves to make Q,an 
0(1) quantity for simplicity in presentation. For 
the tank with PMD sketched in Fig. 4, for 6=0, 
Qfrom eq. (23) is presented for a variety of vane 
lengths V in Fig. 9. The vane edge pinning 
constraint restricts the range of each curve in a 
similar fashion as the curves computed and 
presented in Fig. 5. Q,is maximized for n = 12, V 
= 0.68, which means that the highest rewetting 
flow rate per unit vane length is achieved for 
these conditions for this PMD-type. (It is 
interesting to note that Q,is maximized for n= 12 
and thus a= 15°. This value also corresponds to 
the wedge half-angle yielding the maximum 
capillary flow rate for a fixed volume spreading 
drop 11 .) 


This example optimization is one of several that 
may be constructed for a variety of complex tank 
geometries. Such analytical schemes are quickly 
accomplished, accurate, and trivial in terms of 
commitment compared to numerically based 
techniques. 





Figure 9. Dimensionless flow rate to vane length ratio 
Q.for radial center vane PMD sketched in Fig. 3. 
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Limitations of the Theoretical Approach 


VTRE PMD Rewetting after Thrust Resettlini 


The preceding analysis to compute R assumes a 
priori knowledge of the interior corners of the 
container that satisfy the Concus-Finn condition. 
The analysis also assumes knowledge of a local 
and symmetric equilibrium surface (one of 
perhaps many 10 ). For more complex symmetric 
containers such as those shown in Figs. 4 and 7 it 
is assumed that the interface is also symmetric. 

For the ensuing transient flow problem, the bulk 
interface is assumed to rapidly achieve a constant 
mean radius of curvature R. The interior corners 
must be sufficiently planar such that the flow 
may be approximated by the system defined by 
eq. (1), Fig. 1. The planar interior comers must 
also be of sufficient size such that the interface 
does not pin on wettability boundaries, i.e. the 
terminus of a vane (where the equilibrium 
contact angle is no longer unique). Such pinning 
flows are address analytically by Romero and 
Yost 12 and experimentally by Mann et al. 13 
Slightly non-planar interior 'corners' may be 
treated by a modified analytical approach 14 . 

For cylindrical containers of increasing 
complexity, a generally increasing number of 
constraints must be applied to the solution for R. 
These constraints limit the range of applicability 
of the present solution procedure. Modified or 
alternate techniques may be developed for 
constraint conditions such as edge pinning or 
single interfaces spanning more than one interior 
corner. The more general though complex 
approach of Finn and Neel 9 may also be applied. 
Such techniques will be discussed in a 
subsequent publication as will be the significant 
impact of contact angle hysteresis for real 
systems where d> 0, which has been ignored. 


plication to Tank PMD Rewetting 


The analysis outlined herein naturally applies to 
spontaneous capillary driven flow as occurs in 
liquid propellant tanks following termination of 
thruster firing for orbital maneuvering, docking, 
or tank resettling. Other examples include myriad 
low-g fluids management applications (i.e. on- 
orbit container filling) and drop tower tests. 
Attention here is focused on the former where the 
results of the Vented Tank Resupply Experiment 
(VTRE) provide in-flight data of PMD rewetting 
following thrust resettling. 


VTRE was conducted aboard the Space Shuttle 
in 1996 15 . The experiment explored a variety of 
practical issues concerning propellant 
management in a space-based system. One of the 
tests performed involved thrust resettling of a 
20% filled spherical tank with PMD: 12 axial 
radial (center post) vanes and 12 axial radial wall 
vanes. The test was conducted by exploiting the 
Orbiter primary Reaction Control (RCS) jets to 
settle the liquid contents in a most unfavorable 
location within the tank to observe the 
spontaneous redistribution of the liquid upon 
termination of the thrust. A schematic of the r = 
0.178m tank is provided in Fig. 10a with a cross- 
section in Fig. 10b. The test fluid was R-113 at 
20 °C with <7= 0.0167N/m, p = 7.21- 10' 4 kg/m-s, 
p = 1570 kg/m 3 , and 8=0. 



Figure 10. Spherical VTRE tank with 12 inner and 
outer radial vane PMD: r = 0.1778 m. Tank vent at 
top. propellant outlet at bottom. 



Figure 1 1 . Approximate VTRE interface 
configurations for ~ 20% fill: a. effective equilibrium 
with g > l(\0)' 4 g o (g 0 = 9.8m/s 2 ) acting positive- 
upward, t = 0, b. L(t) during PMD rewetting 
with g-O (10 ~ 6 g 0 ), c. L(t ) at data termination, g-O 
( 10 _6 g o ), t = t f , d. equilibrium, g~ O (10 ~ 6 g 0 ). 
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The equilibrium interface for g ~ 0 is shown in 
Fig. lid — liquid centered over propellant outlet, 
vapor centered over tank vent. During 
unfavorable thruster firing the liquid contents 
reorients to the configuration sketched in Fig. 
11a. Following termination of the thruster firing 
the fluid spontaneously returns to the low-g 
equilibrium configuration of Fig. lid by the 
combined influence of surface tension, surface 
wettability, and container/vane geometry. It is of 
critical design importance to understand 
quantitatively what minimal PMD will produce 
the desired performance. 

As a first application of the theoretical technique 
to model PMD rewetting following termination 
of thruster firing, the VTRE data was re-analyzed 
to determine the transient meniscus tip location 
Lj(t ) in the interior corners of the tank formed by 
the vanes of the PMD. The fact that the tank was 
filled to approximately 20% led to the initial 
condition of a predominantly flat surface (Fig. 
11a) that did not contact the center radial vane 
structure. Thus, upon termination of the thrust, 
rewetting of the tank consisted first of 
spontaneous corner flows along the radial wall 
vanes to the base of the central radial vane 
structure at the propellant exit port. Figs. 1 lb and 
11c. The central radial vane structure was then 
wetted from below and the spontaneous flow 
along this path eventually returned the liquid to 
the equilibrium configuration shown in Fig. 1 Id. 



0.8 1 1.2 1.4 1.6 

sqrt(f) 

Figure 12. Comer tip location L vs. t m for VTRE 
during spontaneous rewetting of two outer wall vanes. 
Solid line is eq. (4) with 5 (vtre = 0.412. 

The flows of interest are identified schematically 
in Figs, lib and 11c. VTRE data for the 
meniscus tip location L of two radial wall vanes 
is provided in Fig. 12. Significant optical 
distortions prevented accurate measurements for 
other vanes. These data are collected from the 


VTRE flight video tapes following termination of 
the Orbiter RCS firing. The Tracker Image 
Analysis System developed by NASA 16 is used to 
digitize the video images. The meniscus location 
is computed by applying optical corrections for 
camera rotation, depth of field, and projection of 
the 3-D spherical flow onto the 2-D CCD array. 
Measurement accuracy is estimated to be better 
than ±5%, the largest uncertainty arising from a 
5% change in scale factor from the front to 
midplane of the spherical tank. A tank flange 
obscured data for time less than that shown on 
Fig. 12. 

For the two vanes analyzed, L(t ) is presented 
against t m in Fig. 12 as suggested by theory, 
eq. (4). The flows are nearly identical reflecting 
the degree of symmetry of the initial condition 
(thrust well-aligned with tank axis) and computed 
slopes for each vane agree to within 4%. Because 
the time for the initial wall rewetting was so short 
(<1.7s), L vs. t also appears linear for this test. 
Nonetheless, the precision of the linear fit for 
L vs. l 1/2 argues favorably for application of 
the transient analysis outlined herein. Thus, 
applying the form suggested by eq. (4) to the data 
of Fig. 12 

W = 0.232 r 1/2 , (24) 

where the experimentally determined coefficient 
0.232m/s 1/2 is accurate to ±5%. Increased 
uncertainty is expected for t < Is. It is insightful 
to mention that for this 0.356m diameter tank 
average corner flow velocities are as high as 
0.232m/s within Is of thrust termination. Such 
velocities increase with container size to the 
Vi-power. Initial velocities in a similar lm 
spherical tank and fluid are likely to be 0.39m/s. 

Substituting the thermophysical properties of 
R-113, eq. (4) is equated to eq. (24) and solved 
to determine %tre = 0.412. This is the 
experimentally determined value of % which 
when used to predict meniscus tip location L(t) 
during rewetting provides the collapse of the 
experimental data illustrated in Fig. 12 and 
prediction by eq. (4) to within ±5%. 

Generalized VTRE Model Section 

Because flight data of PMD rewetting is 
extremely rare it is of value to apply the 
analytical approach of this paper to the VTRE 


NASA/CR— 2002-21 1974 


9 



tank PMD rewetting test despite the fact that the 
spherical VTRE tank with PMD violates 
numerous assumptions: 


1/2 


( 25 ) 


1. The tank is spherical, not cylindrical , and 3-D 
curvature affects might be expected to be 
significant. 

2. The widths of both central and wall radial 
vanes vary with axial location. 

3. The mean radius of curvature R for equilibrium 
interfaces is a significant function of fill level. 

4. The VTRE tank might be considered ‘large’ 
and the rapid formation of a bulk interface with 
constant R seems unlikely. 

5. VTRE experimental data show that the 
rewetting flows along the corners formed by 
the radial wall (outer) vanes eventually pin on 
the vane edges and that single interfaces are 
observed to span two interior corners formed 
by the outer wall vanes. It is noted that both 
occur near the end of the rewetting event. 



Figure 13. VTRE cylindrical tank model: a. cross- 
section identifying inner/outer wetted vanes, b. 
symmetric element of shaded region in a. with JEj 
identifying symmetry planes. 

In the face of such seeming complications the 
analytical technique is applied to model the 
VTRE PMD rewetting event. A generic 
cylindrical model of VTRE is sketched in Fig. 
13a. The smallest symmetrical element is 
depicted in Fig. 13b. Due to the large number of 
vanes, the curved portion of the section is 
approximated as a straight section in Fig. 13b, 
in a fashion after the tank model of Fig. 7 with 
n = 12. However, unlike the example of Fig. 7, 
the curvature of the tank is modeled with 
improved precision by approximating the 
smallest symmetrical element as an isosceles 
triangle, rather than a right triangle. This decision 
is mute for large n, since the difference in models 
is measured by 1 In. Eq. (17) for the cross-section 
in Fig. 13b yields 



and subscripts 1 and 2 denote inner and outer 
vanes, respectively. (Note that a\ = /till.) 

The presumed interfacial configuration of Fig. 13 
leads to eq. (25) for the prediction of Hi 
However, other more preferred configurations 
may arise, several of which are anticipated as 
sketched in Fig. 14. One approach to determine 
the transient flow problem for each configuration 
is to first assume the configuration, compute Hi 
for that configuration using eq. (16), and apply 
the transient solutions of eqs. (4)-(6). The surface 
energy of a given interface configuration will 
help identify preferred states, but mathematical 
proof is required to establish if a given 
configuration is indeed unique 9 . 

Concerning the configurations of Fig. 14: Fig. 
14a is the case under consideration. Cases 14b 
and 14c are the limiting cases of interface 
pinning on % and T), respectively. Case 14d is 
the limiting case of a single interface (/ 2 ) wetting 
two adjacent out vanes, Tfy Case 14e is the 
limiting condition of T) intersecting the interface 
(/,) in Tj. The cases of 14f, 14g, and 14h are 
actually different configurations and not limiting 
cases of the sought configuration 14a. Case 14f is 
the condition where I\ wets both Tj and T) and 
cases 14g and 14h arise when a third interface / 3 
is present: 14g when l\ only wets Tj and 14h 
when I\ wets both Tj and Other 

configurations might be considered. For brevity 
in the following discussion, the cases of 14f, 14g, 
and 14h will not be consider despite being 
increasingly probable as Tj and T) approach 1. 
The notation I t , I 2 , and / 3 is used to identify 
interfaces in the inner and outer vanes and 
between Tj and the outer wall, respectively, as 
indicated in Fig. 14. 
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Fig. 14. Some possible interface configurations for the 
symmetric subsection of the VTRE model of Fig. 13. 


The dimensionless mean radius of curvature 
( l/ 2 ) from eq. (25) with 6 = 0 is presented in 
Fig. 15 for the range of possible T) values 
identified on the figure. For the interface 
configuration depicted in Fig. 13, the possible 
values for are at least constrained by: 

1. < 1, V\ may not contact tank wall. 

2. V 2 < 1, V 2 may not contact center post. 
B.^S’FjSineq/coseq , h does not pin on V\. 

4. < V j tan(/z /4 - a, /4) , h does not pin on l/ 2 - 

5. ^ < f x (i_<f 2 ), Vi does not touch/i. 

6. ^ <2sin(ci! 1 /2)tan(^/4-ci! 1 /4), I 2 does not 
span 2 outer vane corners. 

7. ^ < y/2 £ v , ^cannot exceed tank maximum. 

It is important to repeat that the above list is not 
exhaustive. 



Figure 15. %‘V\,‘V 2 ) for tank sketched in Fig. 13. 
Upper horizontal dashed lines imply V 2 pinning, 
diagonal dot-dashed line implies 'l\ pinning, lower 
horizontal dashed lines imply 'V 2 intersects l x . 


VTRE Model Section: Special Case 


For the special case of ‘1A = 0.35, 6=0, eq. (25) 
is solved and presented in Fig. 15 along with 
constraints #3 through #7 identified for this 
VTRE-like cylindrical model. It is observed that 
the limiting constraint is interface pinning on the 
inner vanes (I\ pins on #3) and the curve for 
larger values of % (smaller ‘V x ) is approximate at 
best. Constraints #5 and #6 are coincidentally 
nearly identical for this special case of V 2 and the 
curve for lower values of % is irrelevant since 
the fluid configuration is no longer even closely 
modeled by the schematic in Fig. 13. 



1A 

Figure 16. for model VTRE tank (Fig. 12) with 
V - 2 = 0.35: Constraints are identified by list number for 
this geometry. 
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For the further restricted case of ‘1) = 0.6 and 
‘If = 0.35, the Surface Evolve r 4 algorithm is used 
to compute the full 3-D surface for a cylindrical 
tank of radius r, diameter £>, and cylindrical 
section length L. The cylindrical tank has circular 
disc end caps (lids). The aspect ratio m of the 
cylindrical portion of the tank is defined by LID. 
The computed equilibrium surface is shown in 
Fig. 17 for 6= 8.11° for a tank with aspect ratio 
m = 3, and 52% liquid fill volume. A computed 
cross-section of the smallest symmetrical element 
at the mid-plane of the tank is shown in Fig. 18. 



Figure 17. Surface Evolver solution of VTRE-like 
cylindrical model with Tj = 0.6, If = 0.35, ft = 8.11°. 
Oval voids on the perimeter are dry region of the tank 
wall. 



Figure 18. Smallest symmetrical element cross-section 
of cylindrical VTRE model shown in Fig. 17: a. 
Surface Evolver solution, b. schematic identifying 
parameters. Compare with Fig. 13b. 


Several Surface Evolver- computed values (SE) 
for the container mean radius of curvature If are 
listed in Tables 1 through 4 for comparison with 
values computed using eq. (25). Flolding all other 
parameters fixed. Tables 1-4 list values for Of 
dependent on contact angle 6, vane lengths If 
and If, aspect ratio m, and liquid fill level. 
Nominal uncertainties for the SE results are 
provided. The two techniques to determine Uf are 
in excellent agreement. 

Local STs-computed values for and for the 
respective surfaces adjacent to Tj and If are 
also listed in the tables for each case. These radii 
are computed in the plane bisecting the container 
normal to the cylinder axis (Fig. 18). The 
differences between df for the tank computed by 
eq. (25) and Of, and df computed by SE provide 
a measure of error for the use of eq. (25) arising 
from the infinite container assumption. This error 
might be considered small in light of such low 
aspect ratio m containers. It is clear from Table 2 
that all SE values for Of and Of approach eq. (25) 
values for df as m increases. 


Table 1. Comparison of present theory eq. (25) and 
Surface Evolver (SE) computations: Effect of contact 
angle ft. ‘l\ = 0.6, If = 0.35, m= 1, Qu q = 55%. 


ft 

o 

^theo 

eq.(25) 

%. 

SE 

% 

SE 

err. 

(%) 

SE 

err. 

(%) 

0 

0.1268 

0.1273 

0.138 

7.2 

0.122 

1.6 

5 

0.1271 

0.1276 

0.134 

0.4 

0.122 

1.6 

10 

0.1279 

0.1285 

0.135 

1.5 

0.128 

1.6 

20 

0.1319 

0.1325 

0.140 

4.3 

0.132 

1.5 

30 

0.1399 

0.1406 

0.154 

0.6 

0.141 

1.4 

40 

0.1540 

0.1550 

0.176 

0.5 

0.161 

1.2 

44.7 

0.1639 

0.1653 

0.189 

3.7 

OO 

- 


Table 2. Results of Surface Evolver: Effect of aspect 
ratio m; volume of liquid fixed, ft= 0, ^ heo = 0.1268, 
‘l\ = 0.6, Tj = 0.35. Case m = 2 almost uncovers lid, 
case m = 4 uncovers lid. 


m 

(L/D) 

SE 

% 

SE 

err. 

(%) 

% 

SE 

err. 

(%) 

Q 

(%) 

0.75 

0.1260 

0.158 

3.3 

0.144 

6.9 

73 

1 

0.1273 

0.138 

7.2 

0.122 

1.6 

55 

2 

0.1269 

0.128 

0.0 

0.122 

8.2 

28 

4 

0.1059 

0.107 

0.4 

0.103 

9.7 

14 
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Table 3. Results of Surface Evolver : Effect of liq. fill 


Qu q ; e 

= 0, 

= 0.1268, 

O 

II 

.6, V 2 = 

0.35, m 

Qliq 

%. 

% 

err. 


err. 

(%) 

SE 

SE 

(%) 

SE 

(%) 

30 

0.1273 

0.130 

2.3 

0.122 

8.2 

55 

0.1273 

0.138 

7.2 

0.122 

1.6 

70 

0.1271 

0.144 

2.1 

0.136 

7.3 

80 

0.1260 

0.160 

1.2 

0.153 

6.5 


Table 4. Results of Surface Evolver. Effect of vane 
size/ratio V 2 : @ = 20°, m = 1, Q Uq = 55%, SE errors 
< 2%. 


Ti 

VI 

%.heo 

eq. (25) 

%. 

SE 

% 

SE 

% 

SE 

0.60 

0.35 

0.1268 

0.1325 

0.140 

0.132 

0.60 

0.45 

0.1134 

0.1193 

0.131 

0.125 

0.90 

0.35 

0.1031 

0.1001 

0.103 

0.100 


Comparison of Theory and Experiment: VTRE 

As previously mentioned, the VTRE rewetting 
event only involved the outer radial wall vanes 
due to a low fill level in the spherical tank as 
depicted in Fig. 11. Thus the cylindrical tank 
geometry discussed in this paper that models the 
spherical VTRE tank rewetting event following 
thrust resettling is that of Fig. 7. % for this 
cylindrical model was solved as a function of V 
and presented in Fig. 8. 

By equating radii of the spherical VTRE tank 
and cylindrical VTRE model, and by evaluating 
V based on initial interface location (refer Fig. 
11a) and detailed VTRE design drawings 17 
represented only schematically in Fig. 10, a value 
of V= 0.21 may be determined for the rewetting 
event. As demonstrated in Fig. 8, with n = 12, 
this low value for V shows that, at equilibrium, 
the interface pins on the vane edges and single 
interfaces cover two interior corners formed by 
adjacent vanes. Thus, both constraints #1 and #2 
are violated. Nonetheless, observations of the 
flight video show that such constraints are not 
exceeded during the larger portion of the 
transient event. If these constraints were ignored 
for the transient rewetting one might simply use 
the value of ^computed from eq. (19) with n = 
12 and ‘V = 0.21. As shown using dashed lines in 
Fig. 8, $= 0.396 computed in this manner, which 
is in surprisingly favorable agreement (< 4%) 


with %tre = 0.412 determined experimentally. 
An even better prediction is possible using eq. 
(25) setting ‘fj = 0, Fm\ = 0, with ‘rf = 0.21. For 
this case 3^= 0.419. This value is within < 2% of 
% TRE , the improvement arising from the 
approximation of the symmetric subsection as an 
isosceles triangle as opposed to a right triangle. 
Both predictions, using eqs. (19) or (25), are 
correct to within the experimental uncertainty of 
5% for 4-2 = 0.21. 


Table 5. Predicted and measured ^for VTRE. 


Technique g_ 

Predicted, eq. (19) 0.396 

Predicted, eq. (25) 0.419 

Measured eq. (24), Fig. 12 0.412 


Further Considerations 

Following an acceptable agreement for 
between theoretical predictions and VTRE flight 
results compared in Table 5, the theoretical 
approach, which allows the closed form 
calculation of the most important flow 
characteristics such as rise height and flow rate, 
can be used to compute transient interface shapes 
throughout the container. For example, the 
surface within the smallest symmetric sub-section 
of the cylindrical 24 vane VTRE model (Fig. 13) 
is computed in Fig. 19 at various times. The tip 
rise height and receding bulk meniscus location 
may be determined explicitly by eqs. (4) and (6), 
respectively. The latter is exaggerated by a factor 
of 2 in Fig. 19 to clearly illustrate the ‘draining’ 
of the container by the corner flows. 

The full VTRE model is computed and shown in 
Fig. 20 at time t = 2.5s — the approximate 
duration of the initial VTRE PMD rewetting 
event had all the vanes been wetted. The VTRE 
model with only exterior vanes wetted is also 
computed and shown in Fig. 21 at time t = 2.5s — 
the model of the PMD rewetting process actually 
achieved on-orbit. (Note that t^in Fig 21 without 
the central vanes is significantly larger than ^ in 
Fig. 20 with the central vanes.) Computations of 
such surfaces serve well to illustrate the wealth of 
information contained within the closed form 
analytic solutions reported herein. 
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Concluding Remarks 



Fig. 19 Transient corner flows in smallest symmetrical 
sub-section of 24- vane VTRE model (times t = 0. 0.5, 
1.0, 1.5, 2.5). 



Fig. 20. Full VTRE model surface at t =0. 0.5. 2.5s. 12 
Inner and 12 outer vanes not shown. 


- ~i — j# 



Fig. 21. VTRE model surface without central vanes at 
t = 0. 0.5. 2.5s. Model approximates flight condition 
of PMD rewetting after thrust resettling. 12 outer 
vanes not shown. 
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The literature reports an accurate analytical 
solution approach to predict spontaneous 
capillary flows in containers with interior 
corners. Such flows are important to a variety of 
low-g fluids handling operations including 
propellant management. In this paper a procedure 
is outlined and demonstrated that culminates in 
the prediction of transient flows in complex 
cylindrical containers that are symmetric, or 
where the contact angles 9 t around the projected 
cross-section may be specified. The general steps 
are: 

1. Identify the interior corners of the tank 
satisfying the Concus-Finn wetting 
condition. 

2. Derive the mean radius of interfacial 
curvature ^ for the tank. 

3. Identify and derive the constraints on % 


4. Compute Hj from Hi for each wetting corner 
of the tank and compute important transient 
quantities such as flow distance, flow rate, 
receding meniscus location and entire 
surface shape. 

In this paper the important unknown quantity is 
% the dimensionless mean radius of curvature of 
the interface at equilibrium, knowledge of which 
enables the determination of the correct initial 
condition for the sought transient solutions. The 
theory of de Lazzer et al 3 to compute % is 
modified to account for symmetry planes within 
complex cylindrical tanks. Three cylindrical 
vaned tank-types of increasing complexity are 
modeled to demonstrate the approach to compute 
% a circular tank with central radial vanes 
(Fig. 4), a tank with wall mounted radial vanes 
(Fig. 7), and a combination tank which serves as 
a model for the Vented Tank Resupply 
Experiment (VTRE) Shuttle flight tests (Fig. 13). 
It is shown that even for the most complex tank, 
agreement in Hi for the present theory with 3-D 
numerical predictions is typically better than 5% 
for aspect ratio containers of about 1 or greater. 
The results apply in general to symmetrical 
polygonal tanks and certain tanks with curved 
walls as demonstrated. Some of the limitations of 
the theory are noted. 
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The results of the analysis greatly speed and 
simplify calculations of capillary driven flows in 
complex containers which model important 
problems such as PMD rewetting following 
thrust resettling. Experiments concerning PMD 
rewetting were conducted during VTRE testing 
and these data are digitized and presented in Fig. 
12. Hi computed from the VTRE experimental 
data agrees to ±4% with ^computed using the 
theoretical approach as shown in Table 5, despite 
the apparent violation of a significant number of 
assumptions. 

The all-analytical approach espoused herein may 
be used to quickly and accurately determine 
solutions to problems commonly thought to 
require extensive 3-D transient CFD. The 
approach is ideal for design optimization and an 
example problem is solved. The technique may 
be applied as a guide to CFD modeling, or serve 
as a ‘benchmark’ to numerical techniques in 
certain limiting cases. The analytic approach may 
also be exploited to design test tanks mimicking 
the smallest symmetrical sub-section of larger 
tanks for ground tests (i.e. low-g aircraft). The 
tanks may be significantly smaller than the full- 
scale tanks, or even scale models, making data 
taken from brief periods of low-g more 
representative of on-orbit performance. 
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